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Abstract 

I show that the conventional formulations of lattice domain-wall 
fermion with any finite Ng (in the fifth dimension) do not preserve the 
chiral symmetry optimally and propose a new action which preserves 
the chiral symmetry optimally for any finite Ng. 
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Currently, there are many interesting physics issues in supersymmetry, su- 
pergravity and superstring theory that require nonperturbative (numerical) 
studies, in addition to those longstanding ones in QCD. A viable approach is 
to formulate these theories on a spacetime lattice with domain-wall fermions 
(DWF). The basic idea of DWF jlll2] is to use an infinite set of coupled Dirac 
fermion fields [ips{x),s G (—00, 00) ] with masses behaving like a step function 
m(s) = m9{s) such that Weyl fermion states can arise as zeromodes bound 
to the mass defect at s = 0. However, if one uses a compact set of masses, 
then the boundary conditions of the mass (step) function must lead to the 
occurrence of both left-handed and right-handed chiral fermion fields, i.e., a 
vector-like theory. For lattice QCD with DWF in practice, one can only 
use a finite number (A^^) of lattice Dirac fermion fields to set up the domain 
wall, thus the chiral symmetry of the light fermion field is broken. Now the rel- 
evant question is how to construct the couplings between these A^^ lattice Dirac 
fermion fields such that the exact chiral symmetry can be preserved optimally, 
or in other words, the residual mass of the quark field is the minimal. Simi- 
larly, in numerical studies of A/" = 1 supersymmetric SU (n) Yang-Mills theory 
with Ng domain-wall fermions, it is vital to implement the chiral symmetry of 
the gaugino optimally such that the supersymmetry can most easily emerge as 
an "accidental" continuum symmetry on the lattice. In this paper, I discuss 
how to preserve the chiral symmetry of the light fermion fields optimally, for 
any finite A^s (in the fifth dimension). 

First, we examine the domain- wall fermion action^ with open boundary 
conditions [4 , in the context of lattice QCD, 

Ns 

Adwf = 'Yi'{x,s)[Dy,{x,x')Ss,s' + S^^x'D5{s,s')]ij{x',s') (1) 

s,s'=l x,x' 

where is the 4D Wilson-Dirac operator with a negative parameter —mo 

4 

Dw = V/^ + W -mo, mo G (0, 2) 
tf,{x,x') = -[U^{x)6xi^x+fi - Ul{x')6x,^^-f,] 
W{x,x') = Yni'^^-^:'-^' ~^^^(^)^^','^^+^^~^l(^')^^','■^-t^] 

and 

^5(5, S') = 5s,s' - P-5s',s+l - P+Ss',s-1, P± = ^(1 ± 75). 

"'^In this paper, we suppress the lattice spacings (a and 05), as well as the Dirac and color 
indices, which can be easily restored. Also, it is understood that one can replace D^u with 
its improved lattice Dirac operator, e.g., including the clover-like terms. 
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The boundary conditions are fixed by P+iIj{x, 0) = P-tplx, A^^ + 1) = 0. The 
quark fields couphng to physical hadrons can be constructed from the left and 
right boundary modes 

q{x) = ip{x,l)P+ + ip{x,Ns)P- 

Then the quark propagator in a background gauge field can be evaluated El 
as 



where 



{q{x)q{y)) 

S = 
T = 
H 



1-75^ 

1+75^ 



1-T 

1 + T 
(1 - H)^' 
(1 + H)^^ 



(2) 

(3) 
(4) 
(5) 



In the limit A^^ — > cxd, 5* — > H/y/JP (the sign function of H), then the 
quark propagator Q is chirally symmetric. 

However, for any finite A^^, (0) does not break the chiral symmetry in the 
minimal way. In other words, Q is not the optimal approximation for the 
sign function of H. This can be shown as follows. 

First, we rewrite (jHI) as the partial fraction 



S 



where 



TT _2_ v^n bi 



2n + 1 (odd) 
2n (even) 



(6) 



sec 



TT 



dj = tan 



TT 



Here the symbol r^'^'™'\x) denotes an irreducible rational polynomial of the 
form 



PnX'' ^ H hPo 



+ go 



, {m>n, pi,qi > 0) 



Note that the coefficients 6; and di in (jHl) are independent of (the ratio of 
the maximum and the minimum of) the eigenvalues of H^. As it will become 
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clear later, this feature already rules out the possibility that i?*^" ^'"^(if^) or 
can be the optimal rational approximation of {H'^) 
According to de la Vallee-Poussin's theorem and Chebycheff's theorem [7], 
the necessary and sufficient condition for r^"''"^\x) to be the optimal rational 
polynomial of the inverse square root function < Xmin < a; < Xmax 

is that S{x) = 1 — ^/x r*^"''™^(x) has n + m + 2 alternate change of sign in 
the interval [xmin, Xmax], and attains its maxima and minima (all with equal 
magnitude), say, 

5{x) = -A, +A, • ■ • , (-1)"+-+2a 
at consecutive points (xj, i = 1, ■ ■ ■ ,n + m + 2) 

Now, for Ng = 2n, 5{x) = 1 — xR^'^~^''^\x^) is non-negative for x > 0. Thus 
6{x) does not have any alternate change of sign for any intervals in (0, oo). 
Similarly, for Ng = 2n + 1, 5{x) = 1 — xR^'^'^''\x^) is positive for < x < 1, 
zero at X = 1, and negative for x > 1. Thus 5{x) has only two alternate change 
of sign for any A^^ = 2n + 1. Therefore, according to de la Vallee-Poussin's 
theorem and Chebycheff's theorem, we conclude that both i?'^"'~^'")(x^) and 
i?*^"''"'^(x^) cannot be the optimal rational approximation for (x^)~^/^, and (jH)) 
is not the optimal rational approximation for the sign function of H . In other 
words, for any finite A^^, the domain wall fermion action does not preserve 
the chiral symmetry optimally, which in fact underlies the essential difficulties 
encountered in lattice QCD calculations with domain wall fermions. 

Note that even if one projects out the low-lying eigenmodes of if jH] (or 
just the boundary term of the transfer matrix p]), treats them exactly, and 
transforms H into one with narrower spectrum (i.e., with a smaller value of the 
ratio ^'i^axl '^'min) ^uch that the chiral symmetry of is improved, however, 
in principle, © still does not satisfy the criterion for the optimal rational 
approximation of the sign function of regardless of the spectrum of H . 

The optimal rational approximation for the inverse square root function 
was first obtained by Zolotarev in 1877 ^Oj, using Jacobian elliptic functions. 
A detailed discussion of Zolotarev's result can be found in Akhiezer's two books 
[TIE]- Unfortunately, Zolotarev's optimal rational approximation has been 
overlooked by the numerical algebra community until recent years. 

For lattice QCD with DWF, the relevant problem is how to construct a 
DWF action such that the operator S in the quark propagator (j2)) is equal to 

S = HRz{H^) 

where RziH"^) is the Zolotarev optimal rational approximation for the inverse 
square root of H^. In general, we have two options for Rz, namely. 



R 



(n,n.) 
Z 



do A l + /^Vc2/ 



+ hyc2i-i 
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and 

rn-l 



^(n-l,n)j.^2^| _ nr=l (1 + ^ 74/) 



nr=i(i + /^V4.-i) 

where h"^ = if^/A^j„, A^j„ i^max) is the minimum (maximum) of the eigen- 
values of H"^, and the coefficients do, d'^, q and are expressed in terms of 
eUiptic functions jTT] with arguments depending on n and ^max/^min- 

Now if one could construct a domain wall fermion action such that the 
operator T in (@)) is replaced with 

then one can solve for {ug} such that the operator 5 (jS)) is equal to 
1 - T _ J HR^z'''\h^), Ns = 2n + 1 



S 



(8) 



Note that one does not have the option to put different weights for H in the nu- 
merator and the denominator of ((Tj), since the optimal rational approximation 
of the sign function of H is equal to H times the optimal rational approxi- 
mation of Obviously, the highest degree n one can obtain with Ns 
flavors is only [^]. 

Nevertheless, it seems to be nontrivial to implement the weights {oJs} into 
the DWF action (jT)) such that ^ can be reproduced. 

Instead of working with the domain wall fermion action (0), 1 consider one 
of its variants 12 , which differs from (P) by replacing 5x^x'D5{s, s') with 

D^{x, s; x', s') = 5x^x>5s,s' + [Dyj - I), ^,^,^,1 P^5 s' ,s+i + {D^j - l),j,^,j,,P+5s' ,s-i- (9) 

Then the quark propagator in a background gauge field can be evaluated 
as 

{q{x)q{y)) = (10) 

where is same as Q except substituting H with = 75-D^. Evidently, 
fllOj) does not preserve the chiral symmetry optimally, the argument is same as 
the case of 0- 

In view of (|7j) and (jH)) , now it is clear how to construct the optimal domain 
wall fermion action on the lattice. Explicitly, it reads 

Ns 

vA= ^ '^i!{x,s)[{l + UJsD^)^^^>6s,s' - {'^ - i^sDui)x,x'P-Ss',s+l 
s,s'=l x,x' 

-(1 - UJsD^)^^^>P+6s',s-l\^ix' , s') (11) 
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with weights 



(12) 



where sn(t>s; k') is the Jacobian elhptic function with argument Vs (|T!^ and 
modulus k' = y'l — X^in/ ^max i^max X^in ^^e maximum and the 
minimum of the eigenvalues of H^), and {ujg} are obtained from the roots 



Ms 



Ns) of the equation 



Sz(u) 



1 - '(W) 



,(n,n) 
j(n— l,n) 







[ 1 - ^B}z''"''{u) = , 
It can be shown that the argument Vg in p2|) is 



Ns 



' 1 + 3A 
(TTA)^ 



2n + 1 
2n 



2K' 



(13) 



where 



A 



M 



l^lsn^Pf^ 



(14) 



(15) 



K' is the complete elliptic integral of the first kind with modulus k', and B is 
the elhptic theta function. From (fT^ . it is clear that A~„j, < Us < A~^„ since 
sn2(;)<l. 

The quark propagator in a background gauge field can be derived as 



{q{x)q{y)) 



1 - IbSopt 



1 + l^Sopt 

where Sopt is same as (jH)) except substituting H with Hyj. 

Since the chiral symmetry of (fTBj) is equivalent to iS^pi - 
due to a finite can be measured in terms of the deviation 



(16) 



= max 

V y^o 



opt 



which has a theoretical upper bound fT3], 2(1 — A)/(l + A), where A is defined 
in (fTi|). a function of A^^ and h = A.^„^./A^j„. In practice, with Ns = 32, one 
should have no difficulties to achieve Az < 10~^^ for any gauge configurations 
on a finite lattice (say, 16'^ x 32 at /3 = 6.0). 
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It is simple to incorporate the bare quark mass niq by adding the following 
terms 



ruq 



[V'(a;,l)(l-^iZ^^).,..P+^(x',iV, 



2mo 

Ns){l - UJNsDu,)x,x'P-ij{x\ 1)], 

to the optimal DWF action (fTTjl . and changing the boundary conditions to 

P+tlj{x,0) = --^P+^(x,iV,), P_^{x,N, + l) = --^P_^(x,l) . 
Znio zrrio 

After introducing pseudofermions (Pauli-Villars fields) with fixed mass 
niq = 2mo, one can derive the effective 4D lattice Dirac operator for the 
internal quark loops as 

D{mq)=r{D, + mq){l + rD,)-' , r = (17) 

ZTTIq 

where rDc denotes the inverse of the massless quark propagator p6p which 
becomes chirally symmetric in the limit A^^ — * oo. The exponential locality of 
D fll7j) (for any rriq and Ng) has been asserted for sufficiently smooth gauge 
background [T7] . 

In the massless limit (m^ = 0) and Ng ^ oo, D (fT7|) is exactly equal to the 
overlap Dirac operator ^1 , and satisfies the Ginsparg- Wilson relation 

D75 + 75/^ = 2D^,D . 

This implies that D is topologically-proper (i.e., with the correct index and 
axial anomaly), similar to the case of overlap Dirac operator. For any finite Ns, 
D is exactly equal to the overlap Dirac operator with (if^)~^/^ approximated 
by Zolotarev rational polynomial. 

From (fT7|) . the valence quark propagator coupling to physical hadrons can 
be expressed as 

{D, + mq)~' = r(l - rmq)-'[D-\mq) - 1] , 

where D~^{mq) can be computed via the 5-dimensional lattice Dirac operator 
of optimal DWF. Evidently, the valence quark propagator of optimal DWF is 
exactly equal to that of the overlap with Zolotarev approximation. Preliminary 
numerical results have demonstrated that the quark propagator of optimal 
DWF with Ns = 2n is precisely equaP to that of the overlap with (if^)~^/^ 
approximated by R^^~^''^\h^). 

^The relative error between these two quark propagators is always less than 10~^ for 
stopping criterion 10~^^ in the conjugate gradient loops. 
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DWF action 


(j(mo = 1.0) 


cr(mo = 1.8) 


Shamir (Q) 


4.3 X 10^5 


1.9 X 10"^ 


Borici (H)) & © 


5.4 X 10-^ 


1.0 X 10-^ 


Improved DWF ^ 


2.4 X 10-6 


3.5 X 10"^ 


Optimal DWF (HH) 


8.8 X 10-9 


3.8 X 10-1° 



Table 1: The precision of chiral symmetry of the massless quark propagator 
in a gauge background on the 8^ lattice at /3 = 6.0, for various domain-wall 
fermion actions with A^^ = 16. (Note that in the limit Ng oo, a for all 
DWF actions.) 

A simple way to improve the chiral symmetry of DWF action (0) is to 
replace with UsDw, 

Ns 

^dwf = Yi){x,s)[uJsDy,{x,x')5s,s' + 5x,x'D^{s,s)]i){x\s') (18) 

s,s'=l x,x' 

where {ug} are given in (fT^ . It is easy to see that in the limit as — > 0, both (fTTjl 
and ()18j) give the same quark propagator ()16|1 with optimal chiral symmetry. 
In Table 1, the precision of chiral symmetry of each DWF action (with Ng = 
16) discussed above is measured in terms of a = maxjj \(V~^j!^ + 'y5V~^)ij\, 
where = (qq) is the quenched massless quark propagator (with one of 
its endpoints fixed at origin) in a gauge background generated with Wilson 
gauge action at /5 = 6.0 on the 8^ lattice, and the ranges of the indices are: 
1 < i < 12 X 8^, 1 < j < 12. The eigenvalues of are bounded as \\{Hw)\ G 
[0.1848,6.5348] for mo = 1.0, while |A(if^)| G [0.0946,5.7484] for mo = 1.8. 
The quark propagators are computed by conjugate gradient with stopping 
criterion 10-^^. Evidently, the improved DWF action p8|l preserves the chiral 
symmetry much better than Shamir's action and Borici's variant, and the 
optimal DWF (fTT|) is the best among these DWF actions. Finally, we note 
that (fTH|) can be easily implemented for machineries already geared to (jT]). 

In summary, the problem how to construct a DWF action such that the 
effective 4D lattice Dirac operator can preserve the chiral symmetry optimally 
for any given finite Ng, has been solved in (jlip . It provides a better under- 
standing of exact chiral symmetry on a finite lattice, as well as the optimal 
way to tackle nonperturbative issues in QCD or supersymmetric QFTs. 
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